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When attempting to understand the role of disorder in condensed-matter physics, one faces severe experimental and 
theoretical difficulties and many questions are still open. Two of the most challenging ones, which have been debated 
for decades, concern the effect of disorder on superconductivity and quantum magnetism. Recent progress in ultracold 
atomic gases paves the way towards realization of versatile quantum simulators which will be useful to solve these 
questions. In addition, ultracold gases offer original situations and viewpoints, which open new perspectives to the 
field of disordered systems. 



C/2 

i 



o 

(N 
> 

\o 

O 



o 



X 



Phase coherence and interference effects underlie many 
basic phenomena in mesoscopic physics, for instance 
electronic conductiorP, magnetism^, superfluidity and 
superconductivity^! or the propagation of light and sound waves in 
inhomogeneous medial Both also play central roles in high-precision 
devices such as interferometers, accelerometers and gyroscopes. In this 
respect, an important issue concerns the effects of disorder, i.e. of small 
random impurities, which cannot be completely avoided in real-life 
systems. A priori, one may expect that weak disorder slightly affects 
most physical systems and that averaging over the disorder smoothens 
possible effects. One may also expect that, in quantum systems, the 
spatial extension of wavefunctions leads to even weaker effects, via a 
kind of self-averaging. In fact, these naive ideas turn out to be wrong. 
Disorder often leads to subtle situations in which strong effects survive 
averaging over the disordeP, in particular in the quantum world. One 
of the most celebrated examples is Anderson localizatiorPl (AL). It is 
now understood that AL results from interference of the many paths 
associated to coherent multiple scattering from random impurities, 
yielding wavefunctions with exponentially decaying tails and absence 
of diffusion^. This strongly contrasts with the Drude-Boltzmann theory 
of classical transport, which predicts that incoherent scattering induces 
diffusiorP. 

Anderson localization was first introduced for non-interacting quan- 
tum particles to explain the absence of electronic conduction in cer- 
tain dirty solids^, but remained elusive for matterwaves. It was real- 
ized later that it is actually ubiquitous in wave physics^, paving the 
way for the first observations of AL, using classical waves, e.g. light 
in diffusive media-33 and photonic crystalline microwaves33 and sound 
waveP. In condensed-matter physics, AL is now considered a funda- 
mental phenomenon underlying certain metal-insulator transitions, but 
complete theory of disordered solids should incorporate Coulomb inter- 
action, the underlying crystal structure, interaction with phonons, and 
magnetic effects. Unfortunately, understanding the physics of even the 
simplest models including all ingredients poses severe difficulties and 
many issues are still unsolved or even controversial. The most challeng- 
ing ones concern the interplay of disorder and inter-particle interactions, 
and spin-exchange couplings. 

Surprisingly enough, atomic physics offers new approaches to these 
issues. The field of ultracold atoms has been developing rapidly in 
the past decades, making it possible to produce, probe and manipulate 
Bos d 13 * 14 ! and FermP^THl gases with unprecedented versatility, tunabil- 
ity and measurement possibilities (Box 1). Control in these systems 
is now such that ultracold atoms can realize quantum s imula tor d 17 l 18 ( 
i.e. platforms to investigate various fundamental modeli^E^l Land- 
mark results have already been obtained, e.g. observation of Mott 
insulator d 23 * 25 ! Tonks-GirardeaiPSEIl Berezinskii-Kosterlitz-Thouless^l 
physics, and magnetic-like exchangd 2 ^ 3 ^. Investigation of Bose-Einstein 



condensates (BECs) in disordered potentialsSIS] has also emerged in a 
quest for direct signatures of AL of matter-waves. Joint theoreticai^EU 
and experimental effort PZM^ made it possible and two groups succeeded 
recently in observing one-dimensional Al l 44 ! 45 ! 

Prima facie, the discovery of this 'Holy Grail' might mean the end 
of a quest. On the contrary, it is just a beginning as the two experi- 
ments of refs. l44l|45l open unprecedented paths to pursue many outstand- 
ing challenges in the field of disordered systems. Direct extensions in- 
clude studies of metal-insulator transitions in dimensions larger than one, 
and of the effect of weak interactions on localization, for which many 
questions are debated. For stronger interactions, single-particle local- 
ization is usually destro yed, but new conce pts such as many-body An- 
derson localizatiorPIEIl and Bose glas£Sl£l provide original paradigms, 
which renew our understanding of these issues. Experiments on ultracold 
atoms with controlled disorder and controlled interactions can also be 
extended to other systems where disorder plays important roles. For in- 
stance, combining spin exchange implementa tioiP^ l and disorder opens 
the route towards random field-induced ordei^^EI] and spin glasse^^2l. 
These few examples illustrate all the promises of an emerging field, i.e. 
quantum gases in controlled disorder. In this paper, we review theoretical 
and experimental progress in this area and discuss perspectives that are 
now within our grasp. 

The nature of Anderson localization 

Localization, as introduced by P.W. Anderson in 1958, is strictly speak- 
ing a single-particle effecfS. Consider the wavefunction i//(r) of a free 
particle of mass m and energy E, in a rf-dimensional quenched disor- 
dered potential V(r), which is solution of the Schrodinger equation 
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Eifr(r) = -—— tfr) + V(r)(Kr). 
1m 
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While in free space, \j/{r) is an extended plane wave, it can be shown 
rigorousl y 1 60 ! 61 ! that, in the presence of disorder, any solution with arbi- 
trary E is exponentially localized in ID, i.e. ln(|(/^(z)|) ~ \z\/L\ oc , with 
localization length L loc (E) oc / B , where / B is the transport (Boltzmann) 
mean-free path. Eventhough L loc often increases with E, it is striking 
that interference effects of multiply scattered waves is strong enough to 
profoundly affect iA(z), even for very high energies. In 2D, the situation is 
simila^l but interference effects are weaker, and L loc oc / B exp(7rW B /2) 
where k = ^2mE/ti would be the particle wavevector in free space. 
Hence L\ oc increases exponentially for k > 1// B , inducing a crossover 
from extended to localized states in finite-size systems. The situation 
differs dramatically in 3D where a proper phase transition occurs at the 
so-called mobility edge k mob : While low-energy states with k < k moh are 
exponentially localized, those with k > k mob are extended. The exact 
features of the mobility edge are unknown, but approximately captured 
by the Ioffe-Regel criteriorPSI, which basically states that localization 
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requires the coherence volume contain several scattering processes. In 
other words, coherence must survive on longer distances than the mem- 
ory of the initial particle direction, thus yielding k mob ~ 1// B . 

Direct observation of Anderson localization of matter- 
waves 

Observing AL of matter- waves requires several challenging conditions. 
First, one must use weak-enough disorder so that interference effects at 
the origin of AL dominate over classical trapping in potential minima. 
Second, one must eliminate all perturbations such as time-dependent 
fluctuations of the medium, or inter-particle interactions. Finally, one 
must demonstrate exponential localization, not only suppression of trans- 
port as it can also arise from classical trapping. While these conditions 
are very demanding in condensed-matter physics, they can be accurately 
fulfilled with ultracold atoms, using i) controlled disorder, ii) negligi- 
ble interactions, iii) strong isolation from the environment, and iv) direct 
imaging of atomic density profiles. This way, direct signatures of AL 
of non-interacting matter-waves were reported in refs. |44]|45] As we 
shall see, these two experiments are complementary rather than similar 
because they significantly differ as regards both observation scheme and 
class of disorder. 

In ref. [44] a weakly interacting BEC is created in a trap, which is 
abruptly switched off at time t = 0. Then, the condensate expands in 
a guide and in the presence of disorder (Fig. Q}i), created with optical 
speckle (Box 2). This physics is captured by the Gross-Pitaevskii equa- 
tion 

ih-f = -——iji + V(r)</r + (#|V, (2) 
at 2m 

which corresponds to Box 1 Hamiltonian (1) in meanfield r egime . The 
dynamics of the BEC can be understood in a two-stage schemd^H. First, 
it is dominated by interactions and the BEC expands, creating a coherent 
wavefunction with a stationary momentum distribution, D(k) oc 1 - (fa;) 2 , 
where £ = %j ^jAm^i is the initial healing length, which measures the ini- 
tial interaction strength^. Second, once the expansion has strongly low- 
ered the atomic density |i^(z)| 2 , the interaction term vanishes and we are 
left with a superposition of (almost) non-interacting waves tf/ k , the pop- 
ulation of each is D(k). Then each if/^ eventually localizes by interacting 
with the disordered pot ential , so that ln(|^(z)|) ~ |z|/Li oc (A;), and the total 
BEC density reduces tcPUH n BEC (z) - jj dk D{k)(\tjj k (z)\ 2 ). Direct imag- 
ing of the localized matter-wave reveals exponentially decaying tails®! 
with a localization length equal to that of a non-interacting particle with 
k = 1/f (Inset of Fig. QJ). Hence, this experiment corresponds to a 
'transport scheme', which probes AL of non-interacting particles with a 
wavevector controlled by the initial interaction, via £ . 

In contrast, the experiment of ref. [45] uses to a 'static scheme'. The 
interactions are switched off already in the trap via Feshbach resonances, 
so that the gas is created in a superposition of a few (typically 1 to 3) 
low-energy, single-particle eigenstates. They are subsequently imaged in 
situ, revealing exponentially decaying tails (Fig. |TJ>). It is worth noting 
that ref. |45] uses a ID quasi-periodic, incommensurate lattice (Box 2), 
thus realizing the celebrated Aubry-Andre model^f 

H = - ^ J (a\a, + h.c.) + ^ A cos(2nfij + (/>) a\cij (3) 

0,0 j 

i.e. Box 1 Eq. (2) with U = 0, Vj = Acos(2;r/?/ + cf>), and /? an irrational 
number. Differently from the case of truly disordered potentials, there is 
a metal-insulator transition (mobility edge) in ID, which is theoretically 
expected at A/ J = 2. 

These works open new horizons to further deepen our knowledge of 
AL in various directions. In ID, although all states are localized, subtle 
effects arise in correlated disorder, for instance in speckle potentials^. 
To lowest order in the disorder amplitude, V R = -\/(V(z) 2 ), the Lyapunov 
exponent, y(k) = l/Li 0C (£), can be calculated analytically^ and one finds 
y(k) oc (V(2k)V(-2k)) /k 2 , enlightening the role of coherent second-order 
back-scattering, +k — > - k — » + k, in the localization process. Since the 
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Figure 1 I Experimental observation of Anderson localization of mat- 
terwaves with Bose-Einstein condensates, a) Experiments of Institut 
d'Optique (coutesy of V. Josse and P. Bouyer): An interacting BEC expands 
in a tight 1D guide (in red) in the presence of a speckle potential (in blue). The 
expansion stops in less than 500ms and the density profile of the condensate is 
directly imaged (shown in orange-green; from the data of ref. [44]. The column 
density, plotted in semi-logarithmic scale in the inset, shows a clear exponen- 
tial decay characteristic of Anderson localization. The localization length L loc 
extracted by fitting an exponential exp(-2|z|/L loc ) to t he exp erimental profiles®!, 
shows a good agreement with theoretical calculation s ' 35 ! 36 !, b) Experiments of 
LENS (adapted from ref. 45 with permission of the authors): A non-interacting 
BEC is created in a combination of a harmonic trap and a 1 D bichromatic lattice. 
The plot shows the exponent a of a fit of a function exp(-|(.i - xo)/l\") to the 
tails of the condensate at equilibrium in the combined potential, versus the ratio 
of the disorder strength (A) to the site-to-site tunnelling rate (/). The onset of 
localization corresponds to the crossover to a -> 1 for A/J > 9. The inset shows 
a plot of the density profile of the condensate together with the fit for A/J = 15. 



power spectrum of speckle potentials has a cut-off kc, such that C2(2k) = 
(V(2k)V(-2k)) = for k > k c (Box 2), one finds an abrupt c hange 
(effective mobility edge) in the A'-dependence of y for weak disordei^^ 
While y(k) ~ for k < kc, higher-order scattering processes dominate 
for k > k c and y(k) ~ V£. 

In dimensions higher than one, the self-consistent theory of 
localizatiorPSl allows one to calculate L loc and exhibits a mobility edge 
in 3D. It is however known that it is not fully exact. Therefore, a ma- 
jor challenge for disor dered, ultracold atoms is to extend the works of 
refs.[44][45]to twd 70|7 'l and three^l dimensions. Definitely, observing the 
3D mobility edge would be a landmark result, which may stimulate fur- 
ther theoretical developments and drive new approaches by providing 
precise measurements of the mobility edge k mob and the corresponding 
critical exponents, which are unknown. 
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Box 1 | Ultracold quantum gases. 



Creating and manipulating ultracold gases 

Ultracold quantum gases are dilute atomic systems that are cooled down to temperatures of the order of a few tens of nano- Kelvins and confined 
in immaterial traps using combinations of magnetic and optic fields' 4 16 . Owing to strong dilution, the prominent inter-particle interactions 
are two-body interactions while many-body interactions can often be ignored. At ultra-low temperatures, s-wave scattering dominates and the 
interaction is accurately described by a contact pseudo-potential 13 15 . In the general case of mixtures of atoms in different species (or different 
internal states), the physics is thus governed by the Hamiltonian 



ti 2 V 2 
2m c 



+ VJr) 
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(1) 



where ^ and are the field operator and the mass of an atom of species tr. The first integral in Box 1 Eq. (1) represents the single-particle 
Hamiltonian where the potential Vo-(r) is controlled by the configuration of the magnetic and/or optic fields (Box 1 Fig la). In most cases, it is 
nearly a harmonic trap'41^(Vo-(jc, y, z) = Ti(£{x,y,z) "VW^^ 2 /2), the anisotropy of which can be adjusted in experiments. For instance, making 



it strongly anisotropic offers the possibility to produce one- 



or two- - dimensional gases. Another useful possibility is to create a guide 



for the atoms using a strongly focused laser beam . The second integral in Box 1 Eq. (1) represents the interaction operator where g lT ^i is the 
coupling constant for interacting atoms of same or different species (g atCr > > and g^^ < correspond to repulsive and attractive interactions, 
respectively). Interestingly, the value and the sign of g^ can be controlled in quantum gases using Feshbach resonances 22 . 



Optical lattices 

Considering different limits of Hamiltonian (1) allows one to design various models initially introduced in the context of condensed-matter 
physics, but here in a controlled way. One important example is that of optical lattices, which are produced from the interference pattern of 
several laser beams 20 22 . The matter-light interaction creates a periodic potential whose geometry and amplitude are determined by the laser 
configuration and intensity. Both can be controlled in experiments. For instance, using pairs of counter-propagating laser beams (Box 1 Fig lb), 
the lattice potential reads V^ n (x, y, z) = V° YiM X ,yd cos(2& L f ), where V° is the lattice depth and k L the laser wavevector. In deep lattices, the 
atoms get trapped at the periodically-arranged minima of the lattice potential (so-called lattice sites). They can jump from site to site via 
quantum tunnelling (with a rate J), and two atoms interact only in the same site (with an energy U). This physics is governed by the Hubbard 
Hamiltonian, i.e. the discrete version of Hamiltonian (1): 

H = - Yj J <r [K.j^J + hx -) + v <r,j "IjUo-j + J 2 U ^ h l,K',f^'-i^-i ( 2 ) 

where the sum over (J, I) covers all sites j and their nearest-neighbour sites /, and a^ j is the annihilation operator of an atom in site j. Hence, 
ultracold atoms (bosons or fermions) in optical lattices mimic the Hubbard model, which is widely considered in condensed-matter physics, 
for instance to capture the essential physics of electrons in solids. However, in contrast to condensed-matter systems, Hamiltonian (2) can be 
shown to be exact in the limit of deep lattices, low temperature and low interactions 19 . The parameters Va-j and in Box 1 Eq. (2) can 

be calculated ab initio from the potential Vo-(r) — > Vo-(r) + V^"(r) and the coupling constant in Box 1 Eq. (1) and are thus controllable in 
experiments. 




Box 1 Figure 1 1 Confining ultracold atoms in magnetic and optic traps, a) Harmonic trapping and laser waveguide (coutesy of V. Josse 
and P. Bouyer). Magnetic coils create a nearly harmonic trapping potential at the bottom of which a degenerate quantum gas, surrounded by a 
cloud of thermal atoms, is formed. A focused laser beam which creates an almost one-dimensional waveguide is also represented, b) Optical 
lattice. The interference pattern of pairs of counter- propagating laser beams form a periodic potential (represented here in two dimensions). 
The atoms are trapped in the lattice sites, but they can tunnel from site to site with a tunnelling rate J and interact when placed in the same site 
with an energy U. 
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Interactions versus Anderson localization 

Another outstanding challenge is to understand how interactions affect 
localization, a question that has proved puzzling from the earliest times 
of AlP3 an( j which is still debated. Common belief is that even weak 
interactions destoy localization. Different approaches however provide 
apparently contradicting answers in different transport schemes. For 
instance, recent numerical calculation ;! 74 ! 75 ! suggest that for expanding 
BECs, repulsive interactions destroy AL beyond a given threshold. Con- 
versely, other recent results^predict that localization should persist even 
in the presence of interactions. Finally, in transmission experiments 
(which amount to throw a mono-kinetic wavepacket to a disordered re- 
gion and measure transmission), perturbative calculations and numeri- 
cal results indicate that repulsive interactions decrease the localization 
length before completely destroying localization^. Since a non-linear 
term is naturally present in BECs (see last term of Eq. 10), and can be 
controlled via Feshbach resonances®, transport experiments with inter- 
acting condensates are particularly promising to address this question. 

A different approach to the interplay of interactions and localiza- 
tion consists in considering a Bose gas at equilibrium in a rf-dimensional 
box of volume CI in the presence of interactions and disorder (Fig. O. 
For vanishing interactions and zero temperature, all bosons populate the 
single-particle ground state, \xo). Very weak attractive interactions are 
expected to favor localization by contracting the Bose gas, but also in- 
duce instabilities for moderate interactions, pretty much like for trapped 
BECs®. Conversely, weak to moderate repulsive interactions do not af- 
fect much the stability, but work against localization, by populating an 
increasing number of single-particle states, \x v ). Weak interactions pop- 
ulate significantly only the lowest-energy states. Since they are strongly 
bound in rare, low-energy modulations of the potential, their mutual 
overlap is small. The gas then forms a Fock state, IT) oc rivC^v^'IO), 
where b\ is the creation operator in state h^ v ). The population N v of 
each is determined by the competition between single-particle energy e,, 
and interaction within each state \x r )- This results in the characteristic 
equation of stated, Ng = f de D n (e)(/J - e)P(e), where D n (e) is the 
density of states and P v = 1/ jdr [^ v (r)| 4 is the participation volume of 
\Xy). This state is an insulator with finite compressibility, k = dN/dfi, 
and can thus be refered to as a Bose glas d 49 ! 50 ! It attains particularly 
interesting features in disordered potentials bounded below (i.e. when 
V(r) > V min everywhere), for which Lifshits has shown'' that the rele- 
vant single-particle states are determined by large-scale modulations of 
the potential. Since they are exponentially far apart, the density of state 
is exponentially small, 2)n(e) ~ exp [-c(e - Vmin)"* 2 ]. As one can see, 
the equation of state is determined by both the density of state Dn(e) 
and the localization via P{e) in the Lifshits tail, which leads us to name 
this state the Lifshits-Anderson glasPO. In the opposite limit of strong 
interactions, there are very many populated |y v ), which thus overlap, and 
the above description breaks down. In turn, the gas forms an extended, 
connected (quasi-)BEC of density n(r) = \p. — V(r)]/g, which is well de- 
scribed in meanfield approach^. This state is a superfluid. Finally, the 
intermediate region interpolates between the Lifshits-Anderson glass and 
the BEC regime. Then, the Bose gas separates in fragmented, forming 
a compressible insulator (Bose glass). The characteristic features of the 
fragments can then be estimated in the meanfield frameworlPH. 

The above description is consistent with the idea that even weak 
interactions destroy single-particle localization. In order to gain fur- 
ther insight, it is worth noting that in interacting systems, the relevant 
states are not the single-particle eigenstates, but are of collective nature. 
For interacting BECs, they are Bogolyubov quasi-particles®. One then 
finds that, although the ground state is extended, the Bogolyubov quasi- 
particles are localized 47 ! 48 ! 82 ! . Their localization properties however dif- 
fer from those of Schrodinger particles, owing to a strong screening of 
disorders. i n ID, the Lyapunov exponent of a Bogolyubov quasiparti- 
cle reads T(k) = [S(k)] 2 y(k), where y{k) is the single-particle Lyapunov 
exponent and S(k) = 2(faf) 2 /(l + 2(fc£) 2 ) is the screening function. One 
thus finds that in the phonon regime (k <k the screening is strong 
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Figure 2 | Effect of interactions in disordered Bose and Fermi gases. The 

gas is described using the single-particle (non-interacting) states [*v>. In the 
presence of disorder, these states, which are localized and distributed around 
in a given volume, are represented by the spheres (in red when they are pop- 
ulated). Bose gas: For non-interacting bosons (top, central panel), the ground 
state, \xo), only is populated. Then, attractive interactions (top, left panel) tend 
to contract this state, thus favoring localization. Conversely, repulsive interac- 
tions (top, right panel) work against localization by populating more and more 
\Xv) states. Fermi gas: In the absence of interactions, a gas of N fermions pop- 
ulates the N lowest-energy \x v ) states (bottom, central panel). Then, each state 
tends to extend under the action of attractive interactions as for maximizing the 
overlap between different populated Ly,,) states (bottom, left panel). Conversely, 
for repulsive interactions, they tend to minimize their mutual overlap, then favor- 
ing localization (bottom, right panel). 



and T(k) <K y(k). Conversely, in the free-particle regime (k » l/<?), the 
screening vanishes and T(k) =^ y(k). Hence, surprisingly, localization can 
survive in the presence of strong mean-field interactions. This poses new 
challenges to ultracold atoms: Not only direct observation of many-body 
AL, but also possible consequences on quantum coherence, sound-wave 
propagation or thermalization process. 

Fermi systems and 'dirty' superconductors 

Consider now a Fermi gas, and focus again on the ground state properties 
(Fig.O. In the absence of interactions, the gas of N fermions populates 
the N lowest single-particle states. For low density, short-range interac- 
tions do not play a significant role as the populated states are spatially 
separated. However, for large-enough density, they do overlap. Then, for 
repulsive interactions, each populated state tends to contract to minimize 
its overlap with other populated states, thus favoring localization. Con- 
versely, for attractive interactions, the populated states tend to extend to 
maximize their overlap, thus favoring derealization. Hence strikingly, 
interactions have opposite consequences for fermions and bosons. 

Perhaps even more fascinating is the possibility to study 'dirty' 
Fermi liquids. Experiments with two-component Fermi gases (e.g. 6 Li or 
40 K), with interactions controlled by Feshbach resonances, have already 
significantly advanced our understanding of the so-called BEC-BCS 
crossove i 15 ! 16 !. On the attractive side of the resonance and for weak inter- 
actions, the Fermi superfluid is well described by the Bardeen-Schrieffer- 
Cooper (BCS) theory and formation of spatially extended Cooper pairs 
consisting of two fermions of opposite spins and momenta. On the re- 
pulsive side, pairs of fermions form bosonic molecules, which undergo 
Bose-Einstein condensation. Although disorder should not significantly 
affect pairing, BCS superfl uidity and BEC superfluidity are expected to 
react differently to disordef^UHl. The famous Anderson theorerrPS indi- 
cates that disorder should not affect very much the BCS superfluid ow- 
ing to the long-range and overlapping nature of the Cooper pairs. Con- 
versely, disorder should seriously affect the molecular BEC, enhancing 
phase fluctuations. 
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Box 2 | Creating controlled disordered potentials. 




In atomic gases, disorder can be created in a controlled way. For instance, the so-called speckle potentials are formed as follows, 109 . A 
coherent laser beam is diffracted through a ground-glass plate and focused by a converging lens (Box 2 Fig. la). The ground-glass plate 
transmits the laser light without altering the intensity, but imprinting a random phase profile on the emerging light. Then, the complex 
electric field £(r) on the focal plane results from the coherent superposition of many waves with equally-distributed random phases, and is 
thus a Gaussian random process. In such a light field, atoms with a resonance slightly detuned with respect to the laser light experience a 
disordered potential V(r) which, up to a shift introduced to ensure that the statistical average (V) of V(r) vanishes, is proportional to the 
light intensity, V(r) oc ±(|£(r)| 2 - <|£| 2 », an example of which in shown in Box 2 Fig. lb. Hence, the laws of optics allows us to precisely 
determine all statistical properties of speckle potentials. First, although the electric field £(r) is a complex Gaussian random process, the 
disordered potential V(r) is not Gaussian itself, and its single-point probability distribution is a truncated, exponential decaying function, 
P(V(r)) - e^'lV^r 1 exp(-V(r)/V R )0(V(r)/V R + 1), where ^(V 2 ) = |Vr| is the disorder amplitude and is the Heaviside function. Both 
modulus and sign of V R can be controlled experimentally 31 : The modulus is proportional to the incident laser intensity while the sign is 
determined by the detuning of the laser relative to the atomic resonance (V R is positive for 'blue-detuned' laser light^ESSIS] an d negative 
for 'red-detuned' laser light 37 ! 39 " 2 !). Second, the two-point correlation function of the disordered potential, Ciir) = (V(r)V(O)}, is determined 
by the overall shape of the ground-glass plate but not by the details of its asperities 109 . It is thus also controllable experimentally 11 . There 
is however a fundamental constraint: Since speckle potentials result from interference between light waves of wavelength coming from a 
finite-size aperture of angular width 2a (Box 2 Fig. la) they do not contain Fourier components beyond a value 2k c , where k c = (2n/A L ) sin(a). 
In other words, Ci(2\l) = for |k| > kc- 

Speckle potentials can be used directly to investigate the transport of matter- waves in disordered potentials 37 40 . They can also be su- 
perimposed to deep optical lattices 88 . In the latter case, the physics is described by Box 1 Hamiltonian (2) with V a j a random variable 
whose statistical properties are determined by those of the speckle potential. In particular, V^j is non-symmetric and correlated from site 
to site. Yet another possibility to create disorder in deep optical lattices is to superimpose a shallow optical lattice with an incommensurate 
perioc^ & 3 l 45 l 87 i In this case, Va-j = Acos(2nf3j + <f>), where A and <p are determined by the amplitude and the phase of the second lattice 
and yS = k 2 /ki is the (irrational) ratio of the wavevectors of the two lattices. Although the quantity V a j is deterministic, it mimics disorder 
in finite-size systems 3 3 ' 34 ! 89 ! 90 ! In contrast to speckle potentials, these bichromatic lattices form a pseudo-random potential, which is bounded 
(| Vq-jI £; A) and symmetrically distributed. 




Box 2 Figure 1 | Optical speckle potentials, a) Optical configuration, b) Two-dimensional representation of a speckle potential. 



Strongly-correlated gases in disordered lattices 

Strong interactions are also very important in various disordered sys- 
tems, e.g. superfluids in porous media or 'dirty' superconductors. Metal- 
insulator transitions attain a particularly interesting, but not fully under- 
stood character in lattice systems. In this respect, the Bose-Hubbard 
model, 

UO i j 

is central in condensed-matter physicsPEUfor it forms a tractable model, 
which captures the elementary physics of strongly interacting systems. 
Hamiltonian © describes bosons, in a lattice with inhomogeneous on- 
site energies Vj, which can tunnel between the sites, with rate J, and 
interact when placed in the same site, with interaction energy U. In- 
terestingly, this model contains the most fundamental two phenomena 
underlying metal-insulator transitions. They correspond to the Ander- 
son transitiorP0 in the absence of interactions (V = 0) as discussed 



above, and to the Mott transitional in the absence of disorder (V = 0). In 
systems dominated by repulsive interactions, density fluctuations, which 
are energetically costy, are suppressed, and a Mott insulator (MI) state, 
I^mi) K llj(Sj)"|0), is formed. Then, the number of bosons per site, 
n = \p/U + 1], where [.] represents the integer part, is determined and 
phase coherence between the lattice sites vanishes. Mis are insulating, 
incompressible, and gapped as the first excitation corresponds to trans- 
fer one atom from a given site to another, which costs the finite energy 
U. At the other extreme, when tunneling dominates, the bosons form a 

superfluid state, l^sp) oc \ Yja.j |0), with normal density fluctuations 
and perfect coherence between the lattice sites. This state is gapless and 
compressible. 

In the presence of disorder, a glassy phase is formed, which inter- 
polates between Lifshits-Anderson glass for weak interactions, to Bose 
glass for strong interaction^. The latter can be represented as I^bg) x 
Ylj (al) n i\0) with tij = [(/i - Vj)/U + 1]. This phase is thus insulating but 
compressible and gapless since the ground state is quasi-degenerated, 
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like in many other glassy system 



With the possibility of realizing 



experimentally systems exactly described by Hubbard models (Boxl), 
ultracold atoms in optical lattices offer also here unprecedented opportu- 
nities to investigate this physics in detail, and to directly observe the Bose 
glass, which has not been possible in any system s o far. Two experimen- 
tal groups have made the first steps in this directiorP^HH. The experiment 
of ref. 87 applied a bichromatic, incommensurate lattice to ID Mott in- 
sulators. Increasing disorder, a broadening of Mott resonances was ob- 
served, suggesting vanishing of the gap and transition to an insulating 
state with a flat density of excitations. Intensive theoretical studies have 
been devoted to understand these results, using quantum Monte-Carlo 1 ^! 
and Density Matrix Renormalization Groups techniques. The results of 
ref.[89|suggest that, in the conditions of ref.[87] one should expect a com- 
plex phase diagram with competing regions of gapped, incompressible 
band-insulator, and compressible Bose glass phases. Clearly, novel and 
more precise detection schemes are needed to characterize this kind of 
physics, such as direct measurements of compressibility^ or condensate 
fraction in superfluid, or coexisting superfluid and MI phases. The latter 
has been approached experimentally in ref. [88] where disorder-induced 
suppression of the condensate fraction in a lattice with super-imposed 
speckle was observed. 

One can also investigate the corresponding Fermi counterparts with 
ultracold atoms. These systems are particularly interesting as they would 
mimic superconductors, even better than bosons. In this respect, an out- 
standing challenge is definitely to understand high-Tc superconductors, 
and possibly important effects of disorder in these systems. Consider the 
two-component (<r 6 {f, J,)) Fermi-Hubbard Hamiltonian 



H = - J* + h.c.) + Yj Va-j al ^j (5) 

a;(jJ) <r,j 




Figure 3 | The spin glass problem. An assembly of spins located at the nodes 
of a cubic lattice interacts according to Hamiltonian (6) where the exchange term 
J i,i only is randomly distributed, and can be either ferromagnetic (blue bonds) or 
anti-ferromagnetic (red bonds). The ground state of the system corresponds to 
the spin configuration that minimizes the total energy. The inherent complexity of 
spin glasses results from frustration which appears when the topography of fer- 
romagnetic and anti-ferromagnetic bonds make impossible to fulfil the local con- 
straints all together. In some plaquetes of four sites, local minimization is easy, 
for instance when all bonds are ferromagnetic (left disk) or anti-ferromagnetic 
(central disk). In some others, it leads to frustration, for instance for odd numbers 
of ferromagnetic and anti-ferromagnetic bonds (right disk). In the latter case, at 
least one spin is frustrated, that is its spin orientation is not unique. Hence, frus- 
tration is at the origin of a manifold of metastable states which corresponds to 
configurations with similar energies. 



For weak interactions, we have a Fermi liquid similar to that discussed 
above. For strong interactions and low temperature, T < U, the Fermi 
gas enters a MI state, pretty much like for bosons, but with a single 
(n = 1) fermion per site (either | or J,). Evidence of vanishing double 
occ upanc y and incompressibility have been reported recently in Fermi 
Ml d 24 l 25 l Then, in the presence of disorder, various phases could be 
searched for, such as Fermi glasses. At even lower temperatures, spin 
exchange starts to play a role, and a transition from paramagnetic to an- 
tiferromagnetic insulator phases is predicted for ~ 4J 2 /U in non- 
disordered systems^. Interestingly, the interplay of interactions and dis- 
order might lead to appearance of novel 'metallic' phases between the 
Fermi glass and the MI. Hence, dynamical mean-field theory 92 at half- 
filling predicts that disorder tends to stabilize paramagnetic and antifer- 
romagnetic metallic phases for weak interactions. For strong interactions 
however, only the paramagnetic Anderson-Mott insulator (for strong dis- 
order) and antiferromagnetic insulator (for weak disorder) phases sur- 
vive. 

Simulating disordered spin systems 

In condensed-matter physics, other important paradigm models where 
disorder induces subtle effects are spin systems, described by the Hamil- 
tonian 



S' I+ J*,S< 



s : , 



)-£h r S;, (6) 



a') 



with either random spin exchange, Jjj, or random magnetic field, h,. 
Ultracold gases can also simulate this class of systems, although not as 
straightforwardly as for Hubbard models. Consider a two-component 
(Bose-Bose or Fermi-Bose) ultracold gas in an optical lattice, as de- 
scribed by Box 1 Hamiltonian (2). In the strongly-correlated regime, the 
couplings between the particles can be understood as exchange-mediated 
interactions between composite (bosonic or fermionic) particles. One 



can then map Box 1 Hamiltonian (2) onto Hamiltonian ((6) with fictitious 
spins encoded in combinations of the annihilation and creation operators 
of the composite particles: S j_ x = (A, + Ap/2, S ;!/ = (Aj - Aj)/2z, and 
S & = 1 /2-A T A, which indeed have commutation relations of spinP, and 
complicated but analytical functional dependence of Jjj and h, on the pa- 
rameters of Box 1 Hamiltonian (2). In the presence of disorder, these pa- 
rameters are randonP^EU and one can reach various limiting cases corre- 
sponding to Fermi glass, quantum spin glass and quantum percolation^!. 

Particularly promising is the possibility of simulating spin glasses^ 
(Fig.0, for which only the exchange term, Jjj, is randomly distributed. 
The phase diagram of (even classical) spin glasses, which is not known 
yet, is an outstanding challenge in condensed-matter physics. The na- 
ture of spin glasses is still debated and there exist competing theories: 
The Parisi replica symmetry breaking^ and the Nelson-Huse droplet 
modePTJ. Ultracold atoms might contribute to the resolution of this is- 
sue, not only on the classical level but also on the quantum level since 
they offer original ways of performing quenched averages. Importantly 
with a view towards testing the replica theory, overlap between two spin 
configurations between two (or more) replicas can be measured directly 
by preparing a pair of 2D lattices with the same realization of disordef^l. 
Quenched averages for systems with binary disorder can also be simu- 
lated by replacing the classical disorder variables by quantum 1/2-spins, 
and preparing them in a superposition stated. 

Yet another fascinating possibility is to simulate various random 
field-induced order (RFIO) phenomena in systems with continuous sym- 
metry, such as BECs or XY-spin m odels with 1/(1) symmetry, or Heisen- 
berg models with S U (2) symmetry ! 54 ! 55 !. A prototype modePS, is the 2D- 
XY version of Hamiltonian ((6) with fixed exchange /y but random field 
h;. In the absence of disorder, symmetry leads to strong fluctuations, 
which suppress long-range order, according to the Mermin-Wagner- 
Hohenberg theorem. Disorder distributed in a symmetric way suppresses 
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ordering even more. Surprisingly however, disorder that breaks symme- 
try might actually favor ordering. This model can be implemented within 
Bose-Bose mixtures^^J where random uniaxial hy can be implemented 
using two internal states of the same atom, coupled via a random Raman 
field, 7ifl(r) , I / i(r) tv I / 2(r) + /i.c. In order to break the continuous symmetry, 
one uses a Raman coupling with constant phase, but random strength. In 
lattice systems, RFIO shows up but is limited by finite-size effects, even 
in very large systems^. In this respect, ultracold atoms offer an alterna- 
tive and fruitful route. Indeed, RFIO turns out to be particularly efficient 
in two (or multi-) component BECs in meanfield regime, where the en- 
ergy functional reads AE m dr n[(ti 2 /4m)(V6>) 2 + ftfl(r) cos(0)L with n 
the atomic density and 8(r) = #i(r) - 6*2(1*), the phase difference between 
the two BECs. This is the continuous counterpart of the 2D-XY model. 
Then, RFIO manifests itself as a fixed 9(r) = ± n/2, and thus allows to 
control the relative phase between the components 55 . This is a striking 
example where ultracold atoms can be used not only to simulate classic 
models, but also offer new and fruitful viewpoints to fundamental issues. 

Further directions 

As the reader has probably noticed, we both are very enthusiastic about 
the future development of the field of disordered quantum gases, and 
probably would like that any interesting direction can be pursued. Lim- 
ited size of the present review has not allowed us to discuss them all, but 
let us briefly mention another. 

Two-component (Bose-Bose or Bose-Fermi) mixtures offer an alter- 
native method to create disorder in optical lattices, namely by quenching 
one component in ra ndom sites, so as to form a background of randomly- 
distributed impuritie d 94 * 96 ! Theoretical analysis using Gutzwiller method 
confirms the appearance of incompressible MI and partially compress- 
ible Bose glass phased. The idea of freezing the motion of the second 
species to form random impurities (i.e. classical disorder) can be gener- 
alized to freezing of any quantum stattfSJ. In this case the system does not 
involve any classical disorder, but nevertheless localization occurs owing 
to quantum fluctuations in the frozen state of the second species. 

One can even relax the freezing condition and consider say two 
bosonic species, one of which tunnels much slower than the other, form- 
ing a quasi-static disorder. In a large region of parameters (for repulsive 
inter-species forces), the ground state corresponds to full phase segrega- 
tion. In practice it is marked by a large number of metastable states in 
which microscopic phase separation occurs, reminiscent of emulsions in 
immiscible fluids^. Such quantum emulsions are predicted to have very 
similar properties to the Bose glass phase, i.e. compressibility and ab- 
sence of superfluidity. Such quasi-static or even time-dependent disorder 
effects have been suggested to underlie the quite large shift of the SF-MI 
transition in Bose-Ferm l 100 ! 101 ! and Bose-Bose^SU mixtures. This issue 
was quite controversial and the most recent work suggests that, while 
indeed the fermions tend to localize the bosons for attractive boson- 
fermion interactions, higher Bloch bands play a significant rol d 103 * 105 j 
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